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O ■ Abstract 

m _ 

We show that decaying axino dark matter with R-parity violation can explain 
the observed excess of the ISOGeV gamma-ray line from the Galactic center in the 
^ ■ Fermi data. The branching fraction of the axino decay into monochromatic photons 

$H ■ can be 0(1), and constraints from continuum gamma-rays and the anti-proton flux 

are ameliorated. The Peccei-Quinn scale of 0(10^^ — 10^^) GeV and the R-parity 
violation parameter of O(10~^^ — 10~^^) are cosmologically favored. 



1 Introduction 



Recently, there is increasing evidence of the excess of the 130 GeV gamma-ray hne from 
the Galactic Center (GC) in the four-year Fermi data P[2l[3llll|5l[6l[71|8l|9]. This may be 
interpreted as a signal of the dark matter (DM), which annihilates or decays around the 
GC. An obstacle to construct a model of the annihilating/decaying DM which explains 
the observed gamma-ray line excess is that the branching ratio of the monochromatic 
photon production must be fairly large. It should be larger than around 0.01 [lU l \TI \ \T2\. 
Otherwise, continuum gamma-rays would hide the line gamma, and anti-protons may be 
overproduced. For instance, if the DM annihilation into photons takes place through 
loops of the standard model (SM) particles, it is difficult to realize such a large branching 
ratio [T3]. 

In this letter, we propose a model of the decaying DM which naturally explains the 
gamma-ray line excess without producing too much continuum gammas and anti-protons. 
A supersymmetric (SUSY) axion model [H] is considered to solve the strong CP problem 
in the framework of the minimal SUSY SM (MSSM). The axino, which is a fermionic 
superpartner of the axion, is a suitable candidate of the DM, if it is the lightest SUSY 
particle (LSP). By introducing small R-parity violations, the axino decays into a photon 
plus a neutrino, and the Fermi gamma-ray line excess can be explained. It is stressed that 
the branching fraction of the axino decay into monochromatic photons typically becomes 
O(10)%, and the constraints from the overproductions of the continuum gamma-ray and 
the antiproton are satisfied. This is in contrast to the decaying gravitino DM scenario, 
where the branching fraction of the monochromatic photon production is suppressed \1U\ . 
Moreover, the present scenario is cosmologically favored, because the lightest SUSY par- 
ticle of the MSSM (MSSM-LSP), e.g., the lightest neutralino, decays by the R-parity 
violating effects before the big-bang nucleosynthesis (BBN) begins. This avoids the cos- 
mological problem associated with a late decay of the MSSM-LSP when the gravitino is 



lighter than the MSSM-LSP|^ On the other hand, the morphology of the gamma-ray line 
signature from the GC seems to favor the annihilating DM scenario rather than that of 
the decaying DM [TU] . Although relatively large gamma-ray signals are expected from the 
Galactic halo in the decaying DM scenario, no such excesses have been observed. However, 
since there are potentially large uncertainties in the gamma-ray data and the DM density 
profile around the GC, it is premature to specify the DM model by the morphology [TOtlTQ]. 

In the next section, the axino DM model will be introduced, and properties of the 
model will be explained, particularly paying attention to the R-parity violating effects. 
We consider the KSVZ axion models [SO]- It will be shown that the model can explain the 
gamma-ray line excess. In addition, several cosmological aspects will be discussed. The 
last section will be devoted to the conclusion and discussion. 



Light axino DM with R-parity violations was considered in Refs. (TSl [IHl [IT] . Ref. [T5] considered 
the gravitino LSP with the axino heavier than MSSM-LSP. 
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2 Axino dark matter with R-parity violation 



2.1 Axino decay rate with R-parity violation 

Let us first introduce R-parity violations. In this letter, we consider a bilinear type of the 
R-parity violation [21], which is characterized by the superpotential, 

W = fl^L,Hu, (2.1) 

where Li and Hu are chiral superfields of the lepton doublet and the up- type Higgs doublet, 
respectively. The index i = 1,2,3 denotes the generation, and /ij is a parameter with 
a mass dimension. Here and hereafter, summation over i is implicitly promised. By 
redefining Lj and the down-type Higgs superfield as L[ = Li — eiHa and if^ = H^ + eiLi 
with ei = where /i is the higgsino mass parameter appearing in the superpotential 

as W = fiHuHd, the R-parity violating superpotential (12. ip is eliminated. Hereafter, for 
notational simplicity, the primes on the redefined fields are omitted. After the redefinition, 
the SUSY breaking potential becomes 

- £rpv = B.LiHu + m^HLiHl + h.c, (2.2) 

where Li is a scalar component of the superfield Li. The coefficients are Bi ~ —Bti and 
2 ^ (^2 _ -jg^^ where S, and m?^ represent soft SUSY breaking parameters 



intheMSSM, -C^oit = {BHuHd + h.c.)+m'jjjHd\'^+mjjjHu\'^+ml \Li\'^ + - ■ ■ . Due to the 
R-parity violating scalar potential (12. 2p . sneutrinos obtain non-zero vacuum expectation 
values (VEVs) as 

, , m? TT cos B + Bi sin B , ^ 

{y.) = ^ — -V, (2.3) 

where tan/3 = Vu/vd is a ratio of the VEVs of the up- and down-type Higgs fields, v = 
^/vl + v'^ ~ 174 GeV, and m?. is a sneutrino mass. 

Before proceeding to discuss phenomenological aspects, several comments are in order. 
It is possible to introduce the bilinear R-parity violating soft terms, LiHu and LiH^ in 
addition to (12. ip . before the field redefinition. The coefficients in (12. 2 p then have additional 
contributions, but the following analysis will not be affected as far as the R-parity violation 
is parametrized by the the sneutrino VEV fl2.3p . Next, trilinear R-parity violating terms, 
LLE and LQD, are also generated by the field redefinition. They are subdominant and 
will be ignored in the following study, because the terms are multiplied by the Yukawa 
couplings. 

The sneutrino VEVs (12. 3p induce mixings between the SM leptons and the gauginos. 
The SM neutrinos mix with the bino and the neutral wino, and the SM charged leptons mix 
with the charged winos. Hence, the R-parity violating parameters are constrained. The 
neutrinos obtain masses of m^, ~ 9'^{^)'^ I'^b{w)^ where rrij^^y^^ is a bino (wino) mass [22l 
21]. For gaugino masses of O(IOO) GeV, = i^i) jv < 10^"^ is imposed to satisfy the 
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experimental bound on the neutrino masses. Also, the cosmological B — L asymmetry is 
preserved for < 10^'' [251 ISSl [271 128]. Other constraints are known to be weaker (see 
e.g., Ref. [21] )■ As we will see, the size of the R-parity violation favored by the Fermi 
gamma-ray line excess is much smaller as Ki ~ 10~^^. 
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The R-parity violation destabilizes the LSP. In this letter, we consider the axino LSP 
scenario in the KSVZ axion models [20]. The relevant interaction terms of the axino are 



where Cy and Cw are model-dependent coupling constants of order unity, ay(= gy/in) 
is the fine structure constant of U(l)y, aw{= oll^'^) is that of SU(2), fa is the PQ scale, 
a denotes the axino, B{W°') is the bino (wino), and B^yiW^y) is the field strength of the 
U(l)y (SU(2)) gauge boson. The axino LSP is stable as long as the R-parity is conserved, 
whereas it decays via the operators (12.41) with the gaugino mixings with the SM leptons, 
once the R-parity violation is turned on. 

First, let us consider the case of Cw = (see Sec. 12.21 for an explicit realization). The 
first term in (12.41) provides interactions of a — — 7 and a — B — Z. The R-parity violation 
opens a decay of the axino through the B — Vi mixing as a — )■ z/j7 and a — viZ. In the 
limit of — /i| ^ m^, the axino decay rate becomes 



1287r3 /2 1 2m2 



r(5 ^ ^ Thz-.^ \ cos^ Ow 1 , (2.5) 



where is the axino mass, and 6w is the weak mixing angle. Here and hereafter, 
r(a — )■ 'yUi/Zui) denotes a sum of the partial decay rates into Ui and Ui. The factor 
{gY{i'i)/\/2)/m^ in the parenthesis comes from the bino-neutrino mixing, and cos6]y 
from the U(l)y gauge boson-photon mixing. Similarly, we obtain 

Tia ^ Z.., . ( 4^sm^ «M,Vl - 4 Vl - ^ - (2.6) 



1287r3 /2 I 2m| \ ml \ 2m? 2m 



where mz is the mass of the Z boson. For m^ ^ m^, the branching fractions are given by 
Br(a — )■ 7Z/) : Br(a — )■ Zu) ^ cos^^iy : sin^6'iy- From the above results, the axino lifetime 
is estimated as 

... 8 X 10-.eo [J^)' (^)-\ (2,7) 

where the R-parity violating parameter is defined as k = . 

The two-body decay of the axino into a photon contributes to the monochromatic 
gamma signal of the Fermi observation. If the axino mass is around 260 GeV, the photon 



See e.g., Refs. [inilSnilSI] for models to explain such a tiny R-parity violation parameter. 
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produced by a — )■ 71/ has an energy of about 130 GeV. According to Ref. p!0], the observed 
excess of the gamma-ray hne is accounted for by a decaying DM, when its hfetime tdm 
and the branching ratio Br(DM — )■ ju) are in the range of rDM/Br(DM — )■ ju) = (1 — 
3) X 10^* sec and Br(DM — )■ ju) ~ 0.01 — 1. The astrophysical constraints from the diffuse 
gamma-rays [32] and neutrinos [33] are also satisfied for such a parameter region. In the 
present model, the branching ratio is around 0.8 for Cw = 0, while the hadronic branch 
from a — )■ Zui is sufficiently small. Thus, the lifetime and the branching fraction which 
are required to explain the gamma-ray line excess from the GC are realized by the axino 



Next, let us focus on the case of Cw 7^ and neglect the contribution from the first 
term in (12.41) . i.e., Cy = 0. The second term in (12.41) provides interactions of a — W'^ — 7, 
a — W'^ — Z and a — — W^. The decays, a — )■ z/j7, UiZ and ifW^, proceed by these 
interactions with the mixings of W^^ — z/j and W"^ — if. In the limit of — yu| ^ mz, 
the decay rate of the axino into a photon and a neutrino is given by 



where the factor {^g2{yii j j m-^ in the parenthesis is the mixing between the wino and 
the neutrino, and sin6'vi/ the mixing between and the photon. Similarly, the decay 
rate oia^Zvi becomes 

while that of a — )• kW, which is a sum of the rates of a — )> lfW~ and d — )■ l^W^, is 



where mw is the mass of the W boson, and the factor 92(1^1) /^w represents the mixing 
between the charged wino and the lepton. Thus, we obtain Br(a — )■ 71/) : Br(a — )■ Zu) : 
Br(a — )■ Wl) ~ sin^6'vv' : cos"^ 9w '■ 2 for rria ^ mz,mw- This results in the branching 
fraction of a — )■ 7// of around 0.09. 

In the case of Cy ~ Cw, both the first and second terms in Eq. (12. 4p contribute to the 
axino decay. The decay rates in such a generic case are summarized in App. El As shown 
there, the branching fractions are determined by a combination of {mj^/myy){Cw/Cy) for 
\m^y — — fi\ ^ mz- Fig. [T] shows the branching ratios of the axino decay into ju, 

Zv and Wl as a function of mj^jmy^ for Cw = SCy/5 (left) and as a function of 5Cvi//3Cy 
for m^/m^Y = 0.5 (right). One can confirm that the branching ratio of a — )■ 7Z/ becomes 
0.8 for {mj^/m^){Cw/Cy) — )■ 0, while it becomes 0.09 for large {mj^/my^r){Cw/Cy). In 
the intermediate regime, the branching ratio of a — )• 71/ decreases due to an interference 



DM. 




(2.8) 



r(s ^ Wl,) 



1287r3 /2 




(2.10) 
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Figure 1: The branching ratios of a into 71/ (red), Zi^ (green) and Wl (blue) as a function 
of m^/m^y for 3Cy/5 = Cw (left) and as a function of 5Cvi//3Cy for m^/m^^ = 0.5 
(right). We assume ~ 260 GeV and — — /i| ^ m^. 



effect and eventually vanishes at {rrij^ /m^^)(Cw / Cy) = ay /aw [cf- Eq. ( ]A.7p ]. In most 



of the parameter space, however, Br(a — )• 7Z/) > 0(0.01) and hence the model can explain 
the gamma-line without overproducing continuum gamma-ray and antiprotons. 

2.2 A model of SUSY axion 

Here, we briefly describe an explicit model of the SUSY axion. Let us introduce PQ 
superfields, $ and with PQ charges of +1 and —1, respectively. Also, PQ quarks, Q 
and Q, are added, which have fundamental and anti-fundamental representations of the 
SM SU(3), respectively, and both of which have a PQ charge of —1/2. The superpotential 
is given by 

WpQ = XX{^^-V^) + k^QQ + Wo, (2.11) 

where X is a singlet superfield, A and k are coupling constants, and Wq = m3/2Mp is 
a constant term with the gravitino mass m3/2 and the reduced Planck scale Mp. The 
coupling constants are taken to be real and positive. Including the SUSY breaking terms, 
ml and m|, the relevant terms of the scalar potential are 

VpQ = m||$|2 + m||$|2 + A2|<|)$-yY + A^I^P(l'^P + l'5r) + (2Am3/2l/^X + h.c.), (2.12) 

where we have assumed the minimal Kahler potential, for simplicity. The VEVs of the 
PQ scalars are given by ($) ~ ($) ~ V, which is related to the PQ scale fa as fa = 
a/2(($)^ -|- (^)^), and the PQ quarks obtain a mass of kV. The axion, which is a goldstone 
boson associated with the VEVs of the PQ scalars, has an anomaly-induced coupling to 
the gluon via PQ quark loops, because the PQ symmetry is anomalous under the QCD. 
Thus it solves the strong CP problem. The coupling constants, Cy and Cw, depend 
on assignments of the U(l)y and SU(2) charge on the PQ quarks. If they are a singlet 
under SU(2) but have (opposite) U(l)y charges, we obtain Cw = and Cy 7^ 0. If 
they are embedded in the SU(5) representation, both Cy and Cw are nonzero and satisfy 
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Cw = SCy/5. For instance, if the PQ quarks are embedded in 5 + 5 of SU(5), the 
coefficients become 3Cy/5 = Cw = 1- 

In this model, the axino, that is a fermionic component of a hnear combination of $ 
and obtains a mass of rria = A(X) ~ 7713/2, where the VEV of X is derived from (12.121) . 
Several additional effects can make the axino heavier or lighter. There can be other SUSY 
breaking contributions to the tadpole term of X in (12.121) . which change the VEV of X 
and hence the axino mass. Radiative correction from the Q{Q) loops can also modify the 
axino mass [3lj. In this letter, we assume that these effects slightly reduce the axino mass, 
and the axino becomes the LSP. 

2.3 Cosmology 

In this section, we discuss several cosmological constraints on the decaying axino DM 
scenario. 



2.3.1 Lightest neutralino 

Let us assume that the lightest neutralino is mostly composed of the bino and that it is 
the MSSM-LSP. In the presence of the R-parity violation, the bino decays into Zui, W'^lf 
and hui due to the sneutrino VEV. The decay rate of the bino is given by [32] 

2xlO-'Bec(^)-^(^)-\ (2.13) 



p{RPV) 
B 

In order for the bino decay not to disturb the BBN, i.e., for the bino lifetime shorter than 
0.1 sec, we need k > 10~^^ for = 1 TeV. 

The bino also decays into the axino through the R-parity conserving operators ( 12. 4p as 
i? — )■ 07 and B — )■ aZ. If this dominates the bino decays, the produced axinos may exceed 
the observed DM abundance. In order to avoid the axino overproduction, the production 
rate should be much less than r^?^^''. The decay rate of the bino into axinos with photons 
or Z bosons is totally given by |36] 

B 

The axino abundance produced by the bino decay in terms of the density parameter 
fla = Pa/ Pa, whcrc pa is the present energy density of the axino and pcr is the present 
critical energy density, becomes 

p(PQ) 

B ^ ^ B 

where h is the present Hubble parameter in units of lOOkm/s/Mpc, and fi^ is the bino 
abundance after the thermal decoupling evaluated as if the bino were stable. Since the 
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bino abundance is large in general, the axion abundance becomes too large, unless the 
branching ratio of the bino decay into the axino is suppressed, i.e., r^?^^'* ^ 

2.3.2 Axino and axion 

Axinos are produced by scatterings of the gluons and the gluinos from the thermal bath 
at the reheating. The thermally produced axino abundance is given by [371 EH] 

where is the running QCD coupling constant at the Tr scale, and Tr is the reheating 
temperature after the inflation. This is valid as long as Tr is larger than the axino mass. 
Thus, the axino can be a dominant component of the DM for Tr ~ 10^ GeV(/a/10^^ GeV)^. 
The abundance of the axion coherent oscillation is estimated as [IQ] 

where 6a denotes the axion initial misalignment angle. For fa ~ 10^^ — 10^^ GeV, we need 
9a ^ O(O.l) in order for the axion abundance to be lower than the DM abundance. 



2.3.3 Saxion 

The saxion, a, belongs to a fiat direction in the scalar potential fl2.12p . which satisfies 
= V^. It obtains a mass, m„, from the SUSY breaking effect. 

Let us estimate the saxion abundance. The saxion sits around the minimum during 
the inflation, which is slightly displaced from the low-energy true minimum, and begins 
to oscillate around the true minimum when the Hubble parameter decreases to with 
an initial amplitude of cTj. The abundance of saxion coherent oscillation is 

where p„ is the saxion energy density, and s is the entropy density. Here, we have assumed 
that the saxion oscillation starts before the reheating process of the inflation is finished, 
which is the case for Tr < 10^° GeV(mo-/500 GeV)^/^. The saxion dominantly decays into 
the axion pair. The lifetime becomes [H] 



with = 2(($)^ — (<l)^)//^. Note that the saxion also decays into a pair of the gluons 
with a branching fraction of ~ {g^/Air'^)'^ |42|, but it does not affect the BBN as long as 
fa < 10^^ GeV is satisfied for ~ 500 GeV 
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We assume that the saxion is hghter than twice the mass of the axino. Otherwise, the axino LSPs 
are overproduced by the saxion decay, a — >■ aa. This indicates rria < 520 GeV. 



7 



The axions produced by the saxion decay contribute to the extra effective number of 
the neutrino species, AN^s |131I111I121I151IS1I1I1I1H1I19], which is given by 



where 



10.75 



nl/3 



er, \Mp 



_pa_ 

Prad 



T=T^ 



(2.20) 



(2.21) 



,Prad/ T=Ta 

with Tcr denoting the temperature at which the saxion decays and prad being the radiation 
energy density. It is estimated as 



AN, 



eff 



1.2 



fa 



1015 GeV 



R 



WGeV J V500GeV7 



-1 



(2.22) 



The contribution should satisfy AA^ofr ^1- In other words, the recent claims of the 
existence of the extra light species, AA^^efr — 1 [50], can be explained by the non-thermal 
axions from the saxion decay. We assume ^ ~ 1 for simplicity. 



2.3.4 Combined constraints 

The combined constraints on a plane of {k, fa) are shown in Fig. [2l In the top panel, we 
have taken Cw = 0, Cy = 1, rria = 260 GeV, = 1 TeV, rria = 500 GeV and cTj = fa- 
in the light blue band, the Fermi 130 GeV gamma-ray line excess is explained. On the 
right side of the black dashed line, the bino lifetime is shorter than 0.1 sec, and the decay 
has no effects on the BBN. Above the orange dot-dashed line, we have Q^f^h'^ < 0.1, and 
the axino abundance produced by the bino decay is sufficiently small. Here, Q^h'^ = 10 
is taken as a reference value of the bino abundance as inferred from Ref. [51]. Below the 
red horizontal line, AA^^efr < 1 is obtained from the axions produced by the saxion decay. 
In the figure, Tr is set so that Q,^^^^h'^ = 0.1 is realized for each Above the blue dotted 
horizontal line, Tr > and the axino DM is thermally produced. 

In the bottom panel, we have taken Cw = 1 and Cy = and assumed that the 
wino is the MSSM-LSP with its mass of = 1 TeV. Note that since the thermal relic 
wino abundance is small, there is no constraint from the axino overproduction by the 
wino decayl*^ In both cases, it is found that the Fermi 130 GeV gamma-ray line excess is 
accounted for without suffering from the cosmological constraints for k, ^ 0(10"^^ — 10^^^) 
and fa - 10^3 _ iqu 

Here we briefly mention parameter dependences on these constraints. For larger neu- 
tralino mass, the axino lifetime becomes longer and the light-blue band moves to the 
bottom-right. Also, the axino abundance from the neutralino decay fl2.15p becomes larger 
and the orange dot-dashed line moves to the top-right. For smaller saxion mass or larger 

^^Similar conclusions hold for the case of the stau MSSM-LSP. 
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Figure 2: (Top) Constraints on the plane of n and fa for the bino MSSM-LSP with Cw = 0, 
Cy = 1, riia = 260 GeV, = 1 TeV, m„ = 500 GeV and CTj = fa- In the hght blue band, 
the Fermi 130 GeV gamma-ray line excess is explained. On the right side of the black 
dashed line, the bino lifetime is shorter than 0.1 sec. Above the orange dot-dashed line, 
the axino abundance produced by the bino decay satisfies ^^^^h? < 0.1, where VL^h? = 10 
is taken as a reference. For each /„, Tr is set so that the axino becomes the dominant 
component of the DM. Below the red horizontal line, AiVeff < 1 is obtained from the 
axions produced by the saxon decay. Above the blue dotted horizontal line, Tr > rria is 
satisfied to account for the axino DM. (Bottom) Same as the top panel, but for the wino 
MSSM-LSP with C^;^. = 1, Cy = and rriy^, = 1 TeV. 
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initial amplitude of the saxion, the bound from A^^ce fl2.22p becomes stronger and the red 
line moves to the bottom. For general values of Cw and Cy, the branching ratio of the 
axino decay into photon may become smaller, as shown in Fig. [H and hence the light-blue 
band moves to the bottom-right to make the axino lifetime smaller. 

Before closing this section, let us comment on the gravitino. The gravitinos are also 
produced by scatterings of the gluons and the gluinos from the thermal bath at the reheat- 
ing. If the gravitino is lighter than the MSSM-LSP, then it dominantly decays into the 
axino and the axion, and hence there is no BBN constraint |39] . The nonthermal produc- 
tion of axinos by the gravitino decay is negligible, since the abundance of the thermally 
produced gravitino is less than that of the axino. On the other hand, if the gravitino 
is heavier than the MSSM-LSP, its decay affects the BBN. The parameter range corre- 
sponding to Tr ~ O(IO^) GeV (just below the red horizontal line in Fig. [2]) is constrained 
depending on the gravitino mass and the MSSM mass spectrum [5T] . 

3 Conclusion 

We have proposed the decaying axino DM scenario as a model to explain the Fermi 130 GeV 
gamma-ray line excess from the GC. It is based on the SUSY KSVZ axion model with 
the bilinear R-parity violation. The model realizes a fairly large branching fraction of 
the axino decay into a photon plus a neutrino. It was found that the Fermi excess is 
accounted for while satisfying the other cosmological constraints for fa — 10^^ — 10^^ GeV 
and K ~ 0(10-12 _ iq-H). 

Compared to another well-motivated decaying DM, i.e., the decaying gravitino DM, 
the decaying axino DM typically has a larger branching fraction into the monochromatic 
gamma. The gravitino universally couples to the lepton and the Higgs superfields, and 
hence the gravitino's decay into hu, Zv and Wl cannot be suppressed in the presence of 
the bilinear R-parity violation. Thus, the decaying gravitino DM is severely constrained 
by the observation of the antiproton flux [10]. 

Let us touch on the collider phenomenology. The MSSM-LSP is stable in the detectors 
for the R-parity violation of k ~ 10-^^ (see e.g., Ref. [S2])- Thus, when the MSSM-LSP is 
neutral, the SUSY events would be detected by searching for signals with a large missing 
transverse momentum. If the MSSM-LSP is a charged particle, it leaves a charged track, 
which is a characteristic signal in the detector. The sensitivities of the LHC searches are 
the same as the standard SUSY searches. 

As mentioned in the introduction, the morphology of the observed gamma-ray signa- 
ture is still premature to refute the decaying DM scenario. If the uncertainties will be 
understood in future, such an analysis can be used for distinguishing the DM models par- 
ticularly between the decaying and annihilating DM models. If the former scenario will 
become favored, the decaying axino DM model can be an attractive candidate. 
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A Formulae for the axino decay rate 

In this appendix, we provide the general formulae for the axino decay rate in the presence 
of the R-parity violation. The interaction between the axion and the SU(2)xU(l) gauge 
supermultiplets is given by 

L^^^ Id^BAy^^y^^^^^ [d'eAw^w^ + h.c., (a.i) 

4V27r/„ J "^V^nfa J 

where aY(w) stands for the fine structure constant of U(l)y (SU(2)), fa is the PQ scale, A is 
the axion superfield, A = (a+ia) /\^+^/2a9+F^9'^, and Wb (VV^) is the super symmetric 
field strength of U(l) (SU(2)). The coefficients Cy and Cw are model-dependent constants 
of order unity. In terms of the component fields, the axino interactions become 

where d is the axino, S^^, O^ixu) is the gauge boson of U(l)y (SU(2)), and B and W^" 
denote the bino and the wino, respectively. 

Due to the sneutrino VEV induced by the bilinear R-parity violation, the SM leptons 
and the gauginos mix with each other, and the axino LSP, d, can decay into ^u, Zv and 
Wl. The mass matrix of the neutralino and the neutrino becomes 

f -mzswcp nizSy/Sp -gY{'^i)/V2\ 

m^y mzcwcp -mzcws^ g2{i>i)/V2 

—mzSwCjs mzCwCjs —fj, 

mz-sw-Sfs -mzcwsp —/lO 
\-gY{i>i)/V2 g2{i>i)/V2 / 

for jC^ = -|(Vio)^MjvV'° + h.c. with 7/5° = (3,^"^, H^, H^,Uif . Here, the subscript 
i = 1, 2, 3 stands for the generation, Cys = cos /3, sp = sin /3, cw — cos 6w and sw = sin 9w 
with the weak mixing angle 6w- The upper-left 4x4 matrix is identical to the neutralino 
mass matrix of the MSSM, and the neutrino masses are approximated to be zero. On the 
other hand, the mass matrix of the chargino and the lepton is given by 



M: 



N 



(A.3) 



w 



\/2mwS/j 



Mc= I V2mwcp li -Yl{u^) | , (A.4) 
g2{i>i) m;, 
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for Cm = -{tp-fMc4^+ + h.c. with = {W-, H;^ , llf and i)+ = {W+, H^, /+ Here, 
y/ is the lepton Yukawa couphng constant, and the upper-left 2x2 matrix is identical 
to the chargino mass matrix of the MSSM. Neglecting the small SM lepton masses in the 
phase space, one finds the decay rates as 



r(a ^ GL) 
with 



1287r3 /2 



m: 



G 



1 



m 



G 



2mi 2m 



^ J FcLiCy, Cw, m^, m^) (A.5) 



I ayCy COS OwU^^j^ + awCw sin OwU^^ 



w\ 



-ayCy sin 9wU^.^ + awCw cos 9wU^^y^\ 



n/2 r 



hW 



for (G, L) 
for (G, L) 

for (G, L) 



(A.6) 



where mc {G = 7, Z, W) is a mass of the gauge boson, and U^^.j^, U^.y^ and ^ are the 
mixings that diagonalize the mass matrices of Eqs. (IA.3I) and (1A.4I) . 



In the limit of — yu| 



— yu| ^ m^, the mixings between the gauginos and the 



hig gsinos become irrelevant. Assuming that the R-parity violation is small, U^^j^, ^ViW 



and U^l^ in Eq. (1A.6P are approximated as 



gY{T>i) 



ayGy COS 9 



w 



72- 



— awCw sin 9 



w 



m 



B 



m 



w 



ayCy sm fc'H^^= h Ow^w cos9w 



m 



B 



72- 



m 



w 



for (G,L) = (7, z/i). 



for (G,L) = (Z,!/- 



for (G,L) = (l^,/i), 



(A.7) 



where {ui) is the sneutrino VEV induced by the R-parity violation. 

As can be seen from (]A.7I) , the branching fractions are determined by (m^/ m^) {Gw / Gy] 
When either Gw or Gy is sufficiently small, the decay rates are reduced to the expres- 
sions in Sec. 12.11 If the PQ quarks are embedded in a complete vector-like multiplet of 
the GUT, they become Gw = 3Gy/5. Moreover, if the GUT relation is assumed for the 
gaugino masses, they satisfy rriy^/mj^ = 3aw/5ay at the one- loop level of the renormal- 
ization group evolution. Then, the decay rate of the axino to 71/ is found to vanish. This 
cancellation is seen in Fig. [H 
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